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Abstract. For an upstream supersonic flow past a straight-sided cone in K.^ 
whose vertex angle is less than the critical angle, a transonic (supersonic- 
subsonic) shock-front attached to the cone vertex can be formed in the flow. In 
this paper we analyze the stability of transonic shock- fronts in three-dimensional 
steady potential flow past a perturbed cone. We establish that the self-similar 
transonic shock-front solution is conditionally stable in structure with respect 
to the conical perturbation of the cone boundary and the upstream flow in 
appropriate function spaces. In particular, it is proved that the slope of the 
shock-front tends asymptotically to the slope of the unperturbed self-similar 
shock-front downstream at inflnity. 

In order to achieve these results, we flrst formulate the stability problem as 
a free boundary problem and then introduce a coordinate transformation to 
reduce the free boundary problem into a fixed boundary value problem for a 
singular nonlinear elliptic system. We develop an iteration scheme that consists 
of two iteration mappings: one is for an iteration of approximate transonic 
shock-fronts; and the other is for an iteration of the corresponding boundary 
value problems of the singular nonlinear systems for the given approximate 
shock-fronts. To ensure the well-definedness and contraction property of the 
iteration mappings, we develop an approach to establish the well-posedness for 
a corresponding singular linearized elliptic equation, especially the stability 
with respect to the coefficients of the elliptic equation, and to obtain the 
estimates of its solutions reflecting both their singularity at the cone vortex 
and decay at infinity. The approach is to employ key features of the equation, 
introduce appropriate solution spaces, and apply a Fredholm-type theorem to 
establish the existence of solutions by showing the uniqueness in the solution 
spaces. 



1. Introduction 

We study the stability of transonic shock-fronts in three-dimensional steady po- 
tential flow past a perturbed cone. The steady potential equations with cylindrical 
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symmetry with respect to the x-axis can written as 

(dxipu) + dy{pv) + ^ = 0, 
dxV - dyU = 0, 
together with BernouUi's law: 



(1.1) 



2 7—1 



(1.2) 



where Kqo := 5?^^ + 7p:jP2o^ determined by the upstream flow state at infinity, 
i.e., the density and velocity (uqo, 0), and y is the distance of the flow location 
in M'^ to the x-axis. In (|1.2p . we have used the pressure-density relation: 



so that the sound speed c — p''''"^)/^. 

For an upstream supersonic flow past a straight-sided cone, a shock-front is 
formed in the flow. When the vertex angle of the cone is less than the critical angle, 
the shock-front may be self-similar and attached to the cone vertex. There are two 
kinds of admissible shock-fronts depending on the downstream condition at inflnity 
(cf. Courant-Friedrichs [18], Chapter VI): transonic (supersonic-subsonic) shock- 
fronts and supersonic-supersonic shock-fronts. In this paper, we are interested in the 
stability of the transonic shock-front, behind which the flow is completely subsonic 
(see Fig. [1}. More precisely, for fixed upstream density > at infinity, our 
problem is to understand the stability of self-similar transonic shock-front when 
the speed of the upstream flow velocity (ucx3,0) is large, equivalently, when the 
Mach number M^o ■= ^ is large. 



1 ' 



7 > 1, 



(1.3) 



Supersonic 




Figure 1. A self-similar transonic shock in three-dimensional 
steady flow past a straight-sided cone 



By scaling the state variables (u, w, p) — > (m, u, p): 




) 



(1.4) 
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the corresponding sound speed becomes c = the equations in (jl.ip remain 
unchanged for the new variables (u, v, p), and the BernouUi constant becomes hao '■= 
i + Therefore, without loss of generality, we can drop "~" for notational 

convenience hereafter to assume that Uao = 1, the Bernoulli constant is 

^^00:=^ + ^pL"'- (1-5) 

Then we have 

= o-(t-i) or v=c^ = ^— = p'^-i fl 6) 

Under this scaling, the problem reduces to the stability problem for self-similar 
transonic shock-fronts in transonic flow past a perturbed cone, governed by (|l.ip - 
(|1.2[) with the Bernoulli constant (|1.5|) , when the Mach number M^o of the upcoming 
flow is sufhciently large, or equivalently, the density is sufficiently small. 

Conical flow (i.e. cylindrically symmetric flow with respect to an axis, say, the x- 
axis) occurs in many physical situations. For instance, it occurs at the conical nose 
of a projectile facing a supersonic stream of air (cf. |18|V The study of supersonic- 
supersonic shock-fronts was initiated in Gu [22] , Schaeffer [30] , and Li [M] first for 
the wedge case; also see Chen [TT] [T21 [13] , Zhang [311[3S], and Chen-Zhang-Zhu [TU] 
for the recent results. The stability of conical supersonic-supersonic shock-fronts 
has been studied in the recent years in Liu-Lien [26] in the class of BV solutions 
when the cone vertex angle is small, and Chen [13] and Chen-Xin-Yin [f 7 in the 
class of smooth solutions away from the conical shock-front when the perturbed 
cone is sufficiently close to the straight-sided cone. 

The stability of transonic shock-fronts in three-dimensional steady flow past a 
perturbed cone has been a longstanding open problem. Some progress has been 
made for the wedge case in two-dimensional steady flow in Chen-Fang [16j and 
Fang [12] . In particular, in [THl [IH], it was proved that the transonic shock is 
conditionally stable under perturbation of the upstream flow and/or perturbation 
of wedge boundary. Also see [5] [6] [T] [TS] [311 [32l [33] for steady transonic flow in 
multidimensional nozzles. 

For the two-dimensional wedge case, the equations do not involve singular terms 
and the flow past the straight-sided wedge is piecewise constant. However, for the 
three-dimensional conical case, the governing equations have a singularity at the 
cone vertex and the flow past the straight-sided cone is self-similar, but is no longer 
piecewise constant. These cause additional difficulties for the stability problem. In 
this paper, we develop techniques to handle the singular terms in the equations and 
the singularity of the solutions. 

Our main results indicate that the self-similar transonic shock-front is condi- 
tionally stable with respect to the conical perturbation of the cone boundary and 
the upstream flow in appropriate function spaces. That is, it is proved that the 
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transonic shock-front and downstream flow in our solutions are close to the unper- 
turbed self-similar transonic shock-front and downstream flow under the conical 
perturbation, and the slope of the shock-front asymptotically tends to the slope of 
the unperturbed self-similar shock at infinity. 

In order to achieve these results, we first formulate the stability problem as a free 
boundary problem and then introduce a coordinate transformation to reduce the 
free boundary problem into a fixed boundary value problem for a singular nonlinear 
elliptic system. We develop an iteration scheme that consists of two iteration map- 
pings: one is for an iteration of approximate transonic shock-fronts; and the other is 
for an iteration of the corresponding boundary value problems for the singular non- 
linear systems for given approximate shock-fronts. To ensure the well-definedness 
and contraction property of the iteration mappings, it is essential to establish the 
well-posedness for a corresponding singular linearized elliptic equation, especially 
the stability with respect to the coefficients of the equation, and obtain the esti- 
mates of its solutions reflecting their singularity at the cone vertex and decay at 
inflnity. The approach is to employ key features of the equation, introduce appro- 
priate solution spaces, and apply a Fredholm-type theorem in Maz'ya-Plamenevskii 
[28j to establish the existence of solutions by showing the uniqueness in the solution 
spaces. 

The organization of this paper is as follows. In Section 2, we exploit the behavior 
of self-similar transonic shocks and corresponding transonic flows past straight-sided 
cones, governed by (|l.l|) - (|1.2p with Bernoulli constant (|1.5[) . In Section 3, we flrst 
formulate the stability problem as a free boundary problem, then introduce a coor- 
dinate transformation to reduce the free boundary problem into a fixed boundary 
value problem, and finally state the main theorem (Theorem 3.1) of this paper and 
its equivalent theorem (Theorem 3.2). 

In Section 4, we establish the well-posedness for a singular linear elliptic equa- 
tion, which will play an important role for establishing the main theorem. Theorem 
13.11 In Section 5, we develop our iteration scheme for the stability problem, which 
includes two steps: one is an iteration of approximate transonic shock-fronts; and 
the other is the iteration of the corresponding nonlinear boundary value problems 
for given approximate shock-fronts. In Sections 6-7, we prove that the two itera- 
tion mappings in the iteration scheme are both well-defined, contraction mappings, 
based on the well-posedness theory for a singular linear elliptic equation established 
in Section 4. This implies that there exists a unique fixed point of each iteration 
mapping leading to the completion of the proof of the main theorem, Theorem 3.1. 

We remark that all the results for the case 7 > 1 is valid for the isothermal case 
7 = 1 as the limiting case when 7-^-1, which can be checked step by step in the 
proofs. 
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2. SeLF-SIMILAR TRANSONIC SHOCKS AND CORRESPONDING TRANSONIC FLOWS 

PAST STRAIGHT-SIDED CONES 

In this section, we exploit the behavior of self-similar transonic shocks and cor- 
responding transonic flows past straight-sided cones, governed by (|l.ip - (|1.2p with 
Bernoulli constant (jl.5|) . 

Let the turning angle of the velocity field right behind the self-similar shock-front 
S be (pi and set b = tani^i. Then v = bu for the velocity field {u,v) of the fiow 
right across S. Assume that the angle between S and the upcoming velocity field 
(1,0) is uji and set r = cottJi. Then the Rankine-Hugoniot conditions on S arc 

[pu] = t[pv], ~[v] = t[u]. (2-1) 

Using (|2.1|1 and the relation v — bu, we have 

r br b + T r.\ 

+ T + T t(1 — or) 

Substitute (|2.2p into Bernoulli's law with Bernoulli constant p.Sp and use i/ = 
jp-. Then a direct computation yields 



b + T /(7-1)(1 + 2t/6-t^) 



"^^(-'-)-^--i-67 r' 2(V+;/5)"^ ' (2-^) 

For 7 > 1 and 6 > 0, we have 

F(0,0) = 0, 9^i^(0, 0) = 1 7^ 0. 

Then the implicit function theorem implies that, in a neighborhood of (0, 0), r can 
be expressed as a function of v, that is, there exists a positive constant z/q such that 

T ~ t{v) for V e [0, vq\. 

Furthermore, there exist positive constants ai and a2 such that, for any v G [0, vq], 
we have 

aivi-^ < t{i') < a2V~<^^ ■ (2-4) 

By (|2.2p . we conclude 

u = 0(l)j/^ 0, iJ = 0(1);^^ ^ 0, p = 0(l) asi/^0, (2.5) 
where 0{\) depends only on 7 and b. Thus, 

= = 0(l)i^^ ^ asi/^0, (2.6) 

where q = \/ + v"^ is the flow speed and 0(1) depends only on 7 and b. 

We now analyze the flow field between the self-similar shock-front S and the 
straight-sided cone. Let ujo be the vertex angle of the cone and k = cotwo- Since 
the equations and the boundary conditions are invariant under the scaling (x, y) — > 
{ax, ay), a 7^ 0, we seek self-similar solutions (w, u) = {u,v){a), a = x/y, as in 
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[18] . Then the flow field {u,v) between the shock- front S and the cone y = kx is 
determined by the following free boundary value problem: 

dcrV + ad^u ~ 0, 



(i-^)a..-(^ + .(i-J))a.. + . = o, 



foro-e(r,K), (2.7) 



{u,v) = {us,vs) , on (J = T, (2.8) 

u — = 0, on (7 = K, (2-9) 

where or k is unknown and determined together with the solution, r and 
(us,vs', ps) E^re determined by the shock polar and the flow direction b right be- 
hind the shock- front S which are given in ()2.2p , and the density p is determined by 
Bernoulli's law with Bernoulli constant (11.511. 




Figure 2. Apple curve and shock polar for the self-similar solutions 

By [T^ , there exists a vertex angle luq = ujq (b) of the cone and the corresponding 
self-similar solution (uo,fo)(o'), a G [r, fc], between the shock-front and the cone as 
the solution of the free boundary value problem (|2.7P " (|2.9p . We assume that the 
flow between the shock-front and the cone is subsonic, which is the case when Moo 
is large (equivalently, poo is small). In this case, we employ (|2.7p to obtain 




where q ~ V^-^ + is the flow speed. It is easy to verify that 

uo(ct) > 0, vo{a) > 0, 
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and uo{a), qoi'^), and the Mach number A/o(ct) are strictly decreasing, while vo{cr) 
is strictly increasing, with respect to a. Therefore, we have 

v(t) v(k) 1 1 
b = ——- < — — — — = tanwo, 1-C-, < n < -, 

u[t) U[K) K 

max uo(ct) = uo{t), 

(tS[t,k] 

max wo(cr) = wo(k) < ?io(T) tanwo, 

<j£[t,k] 

max go(cr) < Qoir), 

max Mo{a) < Mo(r) < 1. 

[TG[r,K] 

In the next sections, we develop a nonlinear iteration scheme and establish the 
stability of self-similar transonic shocks under perturbation of the upstream super- 
sonic flow and the boundary surface of the straight-sided cone. 

3. Stability Problem and Main Theorem 

In this section we first formulate the stability problem as a free boundary value 
problem, then introduce a coordinate transformation to reduce the free boundary 
problem into a fixed boundary value problem, and finally state the main theorem 
(Theorem 3.1) of this paper and its equivalent theorem (Theorem 3.2). 

3.1. Formulation of the stability problem. The stability problem can be for- 
mulated as the following free boundary problem. 

Problem I: Free boundary problem. Determine the free boundary S ~ {x ~ 
0(y)} and the velocity field {u,v) in the unbounded domain {0(y) < a; < ip^^{y)} 
satisfying the equations: 

pv 

d^{pu) + dy{pv) H = 0, 



dxV — dyU = 0, 



V in{0(2/)<x<^-i(y)}, (3.1) 



the free boundary conditions on S: 

MM + HH = o, (3.2) 

-H = [uW{y), (3.3) 
and the slip boundary condition on the boundary surface of the perturbed cone, 

v-tp'{x)u = onB, (3.4) 

where the density p can be expressed as a function of the velocity (w, v) by Bernoulli's 
law: 

7-1 



P ^ 



p{q) ^ (^fioo - ^^-^q'^y \ (3.5) 



with q ^ \/u^ + and Kqo = (7 ^ ^)k,oo- 
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The equations in (|3.ip can be rewritten in the matrix form: 

A{U)d^U + B{U)dyU + Ciy)U = 0, (3.6) 

where U = {u, u)^ and 

To solve the free boundary problem (Problem I), we introduce the following 
coordinate transformation: 

to fix the free boundary: 

f ^ - 'i]cotuji = a: - (j){y), 
n^,^: , , , (3.7) 

[ 77 - ^tantJo = y- 'P{x). 

Then the free boundary S becomes a fixed boundary Fi = = r/cotwi}, and the 
domain {4>{y) < x < ip^^{y)} becomes a fixed domain 

il = {r/cotoji < C < r^cotwo} — {{r, 9) : luq < < oji} . 



In transformation (|3.7p . as a function of y is unknown and can be also considered 
as a function of 77 in the following way: 

if^iv) 0(2/(77 cot wi, 77)). 

Then the transformation is written as 

— 77 cot LUl = X — Ipivi), 

77 - 4 tant^o ^y- ^\x). 

In the case that 7/^(77) is known, we can obtain the expression of 4'{y) from (|3.8p . In 
fact, substituting ^ = -qcotuji into (|3.8p . we have 



x = ?/)(?7), y = (1 — tanojo cot ^1)77 + ip o ■ip{ri). 

Thus, 

dy , • 

— = 1 — tan ujq cot uJi + (p tp, 
drj 

where ip' = '^^ ^ and ^ = -itijl^ jn our case, ip' and tp should be small pertur- 
da; d77 

dy 

bations to tancjo and cottJi, respectively. Hence, we have — > 0, and 77 can be 

dr/ 

also expressed as a function of y, i.e. 77 = 77(77). Then 4>{y) ~ iAviv)) is what we 
need. Therefore, we consider the transformation with formulation (|3.8p from now 
on. Then we have 

y — f] — (p{x) — ^ tan Wo = {'p{'->^) ~ xtanwo) + tan ('0(7;) — ?7Cota;i). (3.9) 
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A direct calculation indicates that the Jacobian matrix of the transformation is 

9(^,77) 1 I I — if' (x) {cot LU I — tp) COtUJi—lp 



d{x,y) 1 — tana;o(cot wi — -0) I tanujo — ip' {x) 1 
d{x,y) I 1 -0 — cotcji \ 



where 



and 



and 



Civ) 



1/77^ 

,0 ; 



J^iU; V) = (5(7?; - A{U; iP)D2U - B(U; iP)DiU, 



1 — tan ciJo (cot tJi — ipj 

1 — tanti;o(cot — %p) 
i ' 



(3.10) 



(3.11) 



'9(C,?/) \ (y9'(a;) — tancjo I + -p' {x)[ij: - cot uji) j 
Then, under the transformation, system (|3.6p becomes 

A{U)d(_U + B{U)d,^U + C{ri)U ^ .^{U-^i') inVL, (3.12) 



C(jr,i}j) := I ?7 ?;(??; V') 
^0 

(I?i, D2) := {d^ + tancjQ 9,,, (cotwi - 7/1) 9^ + 9^). 
Since 7/1(7;) = 0(?/(r; cot wi, 7;)), we have 

7/) — {d^y cotuji + dr^y)(j) = (1 — tantJo cotwi + Lp' ip)4>' , 

and the boundary condition (|3.3p becomes 

• 1 - tanwocottji [7;](1 - tanwocotwi) 

I ~ ip' {x)(p' [y) [u\ + ip'{x)[v\ 



With these, the free boundary problem (|3.ip — (|3.4p becomes the following fixed 
boundary problem: 

Problem II: Fixed boundary problem. Determine the functions (C/; 7/) = 
(u, V] ip) in the unbounded domain: 

f2 :— {rjcoiuji < ^ < 77cotti;o} ~ {(^>^) ^ ijJo < ^ < ^1] 

satisfying system ()3.12p and the boundary conditions: (|3.2p and (|3.13p on Fi := 
— 7?cotwi} and p.4p on Fq := = 7/cotti;o}- 
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3.2. Weighted spaces for solutions. Based on the analysis of the self-similar 
transonic shock solutions in Section 2 and the behavior of solutions to elliptic 
equations at infinity, it is anticipated that the solutions have singularity at the 
origin and decay at infinity. Thus, we need the following weighted spaces as posed 
spaces to accommodate the features of solutions to our problem. 

Let 1 < g < oo and < loq < uji < 2tt. Let 

^ {x e : < r < oo,LJo < 6* < t^i} 

be an unbounded sector, where (r, 9) are the polar coordinates. Then the boundary 
of the domain ^ consists of two rays: 

Fo := {.X e . 0^ < r < cx)} , Fi = {.x £ . 9 = uji,0 < r < oo} . 



W^j^i^) as subspaces of u e W"^'f{&) 
with the norms: 



For any k £ R, m — 0, 1,---, we define the following weighted Sobolev spaces 

soiuew::fm: 

\Hr.d)\\w--m = ' (3-14) 

where 

V{r,e) = {t,0) := {\nr,9) (3.15) 

is a coordinate transformation from (r, 9) to {t, 9). 

Define the norms for the trace of u on each ray Fj of the boundary of 2> by 

\\u{r,ujj)\\^rn-i/,,,^^^^ = ||e*''u(t,cjj)||^y„_i/,,,^g^ , i = 0, 1. (3.16) 

It is easy to see that there exists a constant K , independent of u, such that 



lull ™ < K ||u||tv™,,(^) 



Define 

and denote by Cj™j(f^) the space of functions with norm ||-||(-;m 

When q> 2 and m > 1, the well-known Sobolev imbedding theorem implies that 
^(fe)'(^) embedded in C^™^~^(^), i.e., there exists a constant K, independent of 
u, such that 

||w|lc--i(®) < K I|w|Im/^';:;<'(@) ■ (3.18) 
For functions of single variable defined in M.^ , we can also define the following 
similar weighted norms: 

lkWliM/-'(K+) = l|e'=*'^(e*)L™,,(R) ' ll"Wllc('^,(R+) = l|e''*"(e*)||c™(K) ' (^'^^^ 
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3.3. Main Theorem. The main theorem of this paper is the foUowing. 

Theorem 3.1 (Main theorem). Let q > 2, 1 < ^ < 2, and b> 0. Let 

{(1,0; Poo); {uo,vo){a);(pi = ycotwi}, 

with (T G [cot tJi , cot Wo] and h = M°(cota)lj ' /'^''"^ ^ transonic shock solution to (|3.ip 
when the upstream flow (1,0; poo) pasi the straight- sided cone with y = fo{x) — 
cctanwo- Then there exist positive constants vq, Eq, M, and Ms (M and Mg are 
independent of vq and ) such that, if the Mach number Mao is sufficiently large 
so that V := pj^^ = l/AI^ < vq, then, for any < e < Eq and e ^ v^^, there 
exists a unique solution {U{^,ri);ip{ri)) to the fixed boundary value problem (j3.12p . 
(|3.2p . (j3.13p . and (j3.4p satisfying "0(0) = and the following estimates: 

\\U-Uo\\^.^.^^^<Me, (3.20) 

\\iJj-cotLJi\\r,<MsE, (3.21) 

with ll'llp^ := ll'llvi/°'''(ri) ^" ll'llc"(ri)' provided that, if the perturbed boundary y = 
(p{x) of the cone satisfies ip{0) — and 

y{x) - tancjo||c2^^(K+) + y'ix) - tancjo||vyi^,</(R^) < e, (3.22) 

and the perturbed upstream flow field satisfies 

where fig := {77C0t(ciJi + Sq) < ^ < 77cot(a;o — Sq)} for some small Sq > 0, 

Since Iljp,ip or Tl^^ip is invertible, we conclude the foUowing equivalent result from 
Theorem 13.11 



Theorem 3.2. Suppose that the assumptions of Theorem 13.11 hold. Then there 
exist positive constants uq, Eq, M, and Ms (M and Ms are independent of and 
Eq) such that, if the Mach number M^o is sufficiently large so that v :~ p1^^ = 
1/M^ < vq, then, for any < e < and e ^ , there exists a unique solution 
(still denoted by) {U{x,y);(j){y)) to the free boundary problem (|3.ip - (|3.4p . provided 
that, if the boundary surface y = y^{x) of the perturbed cone satisfies <p(0) = and 

\\ip'{x) ~ tanwollc,2^^(K+) + y'i^) - tanwollw^i^.<!(R^) < e, (3.24) 
and the perturbed upstream flow field satisfies 

where fig ;= {ycot(u;i +Sq) < x < y cotlujQ — Sq)} for some small > 0. Moreover, 
the solution {U{x,y); (f>{y)) satisfies (f>{0) = and the following estimates: 



< Me, (3.26) 
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110' — cota;i||_g < Mss, 



(3.27) 



where (f>' ^ ^ and \\-\\g := IMI w«-'(a+) + l|-|lco(K+)- 




Figure 3. Stability of transonic shock- front solutions 

Remark 3.1. The existence of the perturbed upstream flow field U~ satisfying ()3.25p 
can be obtained by blowing up the angular point and then following the standard 
argument as in Li-Yu |25| , since the equations are still quasilinear hyperbolic under 
the transformation. 

Remark 3.2. Estimates p.26p and (|3.27p imply that the downstream flow and the 
transonic shock-front are a perturbation of the self-similar transonic shock solution. 
Hence, the self-similar transonic shock-front is conditionally stable with respect to 
the conical perturbation of the boundary surface of the cone and the upstream flow 
in the function spaces with restrictions on the downstream flow field both at the 
corner and at infinity. 

4. Well-posedness for A Singular Linear Elliptic Problem 

In this section, we establish the well-posedness for a singular linear elliptic equa- 
tion, which will play an essential role for establishing the main theorem, Theorem 

Let < ujQ < LUi < ^ and set 

:= { (x, y) e : < r < oo,luq < 6 < uji} , 
To := {{x,y) eR^ : < r < oo,e = ujq} , 
Ti -.^ {{x,y) eR^ : < r < 00,6 = uJi} , 
where (r, 9) are the polar coordinates in the plane. 
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4.1. Neumann problem for a singular second-order elliptic equation. Con- 
sider the following Neumann boundary value problem in f2: 

Lof dxxf + dyyip + = / in Q, 

Boip := dy(p — tan /jjQdxf = go on Fq, (^-l) 

Blip :~ dx^p — cotujidyip = gi on Fi. 

We have the following proposition. 



Proposition 4.1. Let 1 < q < oo. The operator {Lq, Bq, Bi) defined in (14. 1|) 
realizes an isomorphism from W^^^j{rt) to W^^1{il) x (wl^^^^'^''^ . Moreover, 
we have the following estimate for the solution to problem (|4.ip ; 

Mw^-\.{n) < I< ( \\f\\w°,nn) + J2 llffjll<Ti/-'(R+) )' (4.2) 

where the constant K is independent of tp, but depends only on q and o^q ( actually 
cot ljJq). 

To prove this proposition, we employ a criterion identified by Hartman-Wintner 
[23] for the uniqueness of solutions to the Dirichlet boundary value problem for 
systems of second-order differential equations. For self-containedness, we give a 
brief description here; for more details, see 

Lemma 4.1. Consider the following boundary value problem for the system of 
second-order differential equations for x G R"; 

r x" + Ai{t) x' + A2{t) X ^ forte {to,ti), 

{ (4.3) 

[x{to)^xiti)^0, 

where Ai{t) and ^2(0 cire n x n real matrices. Assume that there exists a matrix 
K[t) such that 

N = (X°)' - A§ - (Iai - K^) {\aJ K") > 0, (4.4) 

where K° = ^{K + K^) and A2 = ^(^2 + A^). Then problem ([O]) has only the 
trivial solution x = Q. 

Proof. Taking the inner product on the equations with x and integrating from t^ 
to ti yields 

ti 

(:x' -x' ~x- Aix' - X ■ A2X) dt 0. (4.5) 

:> 

Since (x ■ Kx)' — 2x' • K^x + x ■ K'x, we have 

(x ■ KxYdt = / (2x' • K°x + X ■ K'x) dt. (4.6) 

Then ^ 

j ' (\x'~ {]-Aj ^ K'^)x\^ +x- Lx)dt = Q, 
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where L K' - A2 ~ {\Ai - A"") {^Aj - /v") 
Similarly, we have 



Combining the above two identities, we obtain 

/ ' (\x' - {\aJ - K^)x\^ + X ■ Nx)dt = Q. 
Since N is positive definite, we conclude x = 0. □ 



Proof of Proposition 14. 1[ Rewriting the boundary value problem (j4.ip in the polar 
coordinates {r,9), we have 

Lo(p ~ {rdr)'^(p + dg(p + rdr^p + cot 9 detp = r^f in fi, 
Bqip = detp = rgo on Tq, (4.7) 

Blip ~ dgtp = rgi on Fi. 

Employing the transformation V in (|3.15p . i.e., P(r, 9) ~ {t, 9) = (In r, 9), we con- 
vert the infinite sector ^l into an infinite strip: := {{t,9) : t E M., luq < 9 < oji}. 
Accordingly, the boundary value problem (|4.7p is converted to the following bound- 
ary value problem in 

dtt(p + deeip + dt(p + cot 9 deip = e^*/ in @, 

de(p = e*5o on Sq, (4.8) 

dg(p = e*gi on Si. 

Applying the Fourier transformation .^t^x with respect to t, we obtain a family 
of boundary value problems with complex parameter A: 

ip" + cot 9ip' + (-A^ + iX)ip = e2*/ 9e (cjq, wi), 

= e^o e = L,o, (4.9) 

ip' = e^gi 9 = LUi. 

We now employ a Fredholm-type theorem. Theorem A.l in Appendix, to find 
that, if the homogeneous problems of (|4.9p (i.e. f = go = gi = 0) have only the 
trivial solution tp = for all A with ImA = — 1, then, for any (/, 5oi5i) such that 
e*/ e W^°'«(^) and g^ e VFi-i/«'«(Sj), j = 0, 1, the boundary value problem gj]) 
in the infinite strip has a unique solution ip such that e^*(p € VF^'^(^). Moreover, 
the solution satisfies the estimate: 

||e-V|U2.,(s,) < ^(||eVLo,,(sj) + J2 Il5jl!wi-i/...(s,))' (4-10) 

i=o,i 

where K is independent of (/s, but does depend on loq because of the coefficient cot 9. 
Then the results of Proposition 14.11 follow. Therefore, it suffices to verify that, in 
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the case that / = gi = 32 = 0, the boundary value problems (14. 9|) with complex 
parameter A, ImA = — 1, have only the trivial solution ip = 0. 
When f = go = gi = 0, we set y ~ (p' to have 

y" + cot e y' + {- csc^ 9 - \^ + iX) tp ^ for 6* e {ujo, wi), 

y = for6l = cjo, (4.11) 

y = for 9 ^ LUi. 

Write y ^ yi+ 2/2! and A = /i - i. Then - A^ + Ai = -^^ + 2 + 3/ii, and (|4.1ip can 
be rewritten as the following boundary value problems of second-order differential 
equations with real coefficients: 



yi 




yi 



for 9 e {ujQ, wi), 

on 9 = loqjUJi. 

(4.12) 




Let K{9) = K^{9) . Then 



TV 



ai9,n) 
, ai9,^l)^ 



where 



a{9, n) :=(tan 9)' - (- csc^ g - + 2) - (i cot 6* - k)^ 

=/t^ + csc^ 9-- cot^ 6* > 0, 
4 

which implies that the symmetric matrix N is positive definite and hence satisfies 
the criterion (|4.4p . By Lemma [4T] we obtain that y = 0. That is, (p = const. Then 
the equations in (|4.9|1 yields that <y9 = for any ImA = —1. This completes the 
proof. □ 

Proposition 14.11 can be directly applied to a special boundary value problem of 
first-order partial differential equations. 

4.2. A boundary value problem of a singular first-order elliptic system. 

Consider the boundary value problem for the first-order system: 

Ad^U + BdyU + CU = F in f], (4.13) 

U -ao^go on Tq, (4.14) 

U ■ ai ^ gi on Ti, (4-15) 
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where U = {u, v)^ , do = (— tano^o, 1)^, ai = (1, — cotwi)^, F = (/i, /2)^, and 



'1 0^ 
,0 1, 



B = 



V 

-1 0, 



^0 

y 

,0 0, 



(4.16) 



To solve this problem, we first construct a function <t> e ^(-i) (^) such that 



A$ = /2 in f7, 
$ = on To, 
$ = onTi. 



(4.17) 



By virtue of [5S], there exists a unique $ with the following estimate: 



(4.18) 



where K is independent of <&. 

Let u = u + 9y$ and v = v — d^^- Then the boundary value problem (|4.13p - 
(|4.15p is reduced to 



U ■ ao = go 
U ■ ai = gi 



in ri, 

on To, 
on Ti, 



(4.19) 
(4.20) 
(4.21) 



where F = (/, 0) ' , / = /i , go = .9o, and gi = gi. Since 



9.$ 



<K{LUo)\\d.M(i\e)\ 

w°-i{v{n)) 

=K{ujo) \\d^mw°„:^in) <K{ujo) \\f2\\w°''^(n) > 



W<i'i{V(Q.)) 



we have 



< i^(tjo) X! \\h\\w'l,^(n) = C'(^o) ll^llvv°'^<'(o) 



By the second equation of (|4.19p . there exists a potential function Lp such that 
V</3 = [dx'p, dytp) = U . Then the boundary value problem (|4.19p — (|4.2ip can be 
reformulated as a boundary value problem of a second-order elliptic equation: 



y 



dyip - tanujQdxf = go on Tq, 
_ dx^p — cotujidyip ~ gi on Fi. 
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Now Proposition 14.11 yields that there exists a unique solution e W^^^^{^) with 
the following estimate: 

j=o,i 



i=o,i 



(4.22) 



where K depends only on wq, but is independent of (p, F, and gj,j ~ 0, 1. Thus 



there exists a unique solution U G (VF,gf (fi))^ to problem (|4.13p — (|4.15p with the 



following estimate: 

||C/||wi.,(n) < ( li^ll<-(o) + E ll5.llw-^/-(r,) )' (4-23) 

J=0,1 

where K is independent of U, F, and ^j, j = 0, 1, but depends only on ujq. 

With the argument above, we obtain the following corollary of Proposition 14. II 

Proposition 4.2. Let 1 < q < oo. Let F £ {W°^1{fl))^ and g, G Wl^^^''^''^{Tj), 
j = 0, 1. Then there exists a unique solution U G {W^^^ (Q,))'^ to the boundary value 
problem (j4.13p — (j4.15p . Moreover, the solution satisfies estimate (|4.23p . 

Applying the continuity method, we can extend this result to a small "perturbed" 
boundary value problem for the first-order elliptic system. 

4.3. A small perturbed boundary value problem for a singular first-order 
elliptic system. Consider the following boundary value problem: 

Ad^U + BdyU + CU ^F in fl, (4.24) 

U-ao^ go on Pq, (4.25) 

U-ai^gi on Pi, (4.26) 

where A, B, and C are 2x2 matrix functions defined on il, aj — (aji,Q!j2)^ are 
vector functions defined on Pj, j = 0, 1. Then we have 

Proposition 4.3. There exists a positive constant e, depending only on the constant 
K on the right side of estimate (|4.23p . such that, if the coefficients of problem 
(j4.24p — (j4.26p satisfy the following conditions: 



{A-A,B-]3) +C~C + y ||a, ,p . < e, (4.27) 



then, for any F G {W'^^'^ {Q?))'^ and gj G l^(o) ^^''''^ ^ j)^ i ~ 0? 1; there exists a unique 
solution U G (ly^p^'(O))^ to the boundary value problem (|4.24p — (|4.26p . Moreover, 
the solution U satisfies the following estimate: 



i=o.i 

where K is independent of U , F , and gj,j ~ 0, 1, but depends only on loq. 
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Proof. Denote the boundary value problems ([413]) — (|4?T5l) and ([424]) — (|4?26l) by 
linear bounded operators T and T respectively from X = {W^^^ (fl))'^ to Y = 



X n W^(o) (r^)- By Proposition 1131 T is invertible and T'^ is 

j=0,l 

also a linear bounded operator. 

Let Ts = {1- s)f + sT, s e [0, 1]. By virtue of we have 

||(r-T)C/||y < lie 11(711^, 

where is a constant. Since f U ^ TsU + s{f - T)U , by Proposition Ol 

\\U\\x<K {\\TMY + \\iT-T)U\\Y)<K\\TMY+Ki\\U\\x. 

Choosing e sufficiently small such that Ke < 1, we have 

||C/||^ <X i|T,C/||^, (4.29) 

where K is independent of U and s £ [0, 1], but depends only on ujo- 

Then, applying the continuity method, Proposition 14.31 ^^d the uniform esti- 
mates ()4.29p . we completes the proof. □ 

5. Iteration Scheme 

Our iteration scheme for the stability problem consists of two iteration mappings: 
One is for an iteration of approximate transonic shock-fronts; and the other is 
for an iteration of the corresponding nonlinear boundary value problems for given 
approximate shock-fronts. 

Let q > 2 and 4'o iv) V cot ui . Define 

^r=U: ^■(0) = 0,||V'-V'o||ri <t}, 

^ ^ . (5.1) 

Or =\SU = (Su,Sv)' : ||((5w,(5u)||,„i,,,n, < r L 



\ =|(5L/ = {Su,Svy : \\{Su,Sv)\ 



1 1^,0)' (^2) 

Let Ms and M are positive constants to be determined later. In order to find the 
perturbed shock solution to the fixed boundary value problem (|3.12p . (|3.2[) . (|3.13p . 
and p.4p of the self-similar shock solution {Uq; C/q~; tpo), our strategy is as follows: 
Let Eg > be a small constant to be determined later and < e < Eq- Given 
an approximate boundary ijj G Smsej solve the nonlinear boundary value problem 
p.l2p . p.2p . and p.4p to obtain a perturbed solution Ujf, of Uq. Then we use one 
of the Rankine-Hugoniot conditions, (|3.13p . to update the approximate boundary 
and obtain new "i/)*: 

^ [v] ( 1 — tan uJo cot uji) 

[u] + ip'{x)[v] ' (5.2) 

V'*(0) = 0. 

This defines an iteration mapping: ■ ip '—^ ip*- To prove Theorem 13.11 it 

suffices to verify that there exist positive constants Ms and eo such that ^s is a 
well-defined, contraction mapping in Smse ^or any < e < Eq- 
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Since the initial value problem (|5.2p is easier, we will focus mainly on the non- 
linear boundary value problem ()3.12|) . ()3.2p . and (|3.4p for given ip e '^Mse, which 
requires another nonlinear iteration: For given 6U G Omei '•^ linearized boundary 
value problem will be solved in the weighted Sobolev space Wl^{n) to obt am a 
unique solution 5U<, that is defined as an iteration mapping ^ : SU i— > 5C/,. By 
showing that there exist positive constants M and Eq such that ^ is a well-defined 
contraction mapping in Oa/e for any < £ < Eq; we conclude that the nonlinear 
problem (|3.12p . (|3.2p . and l|3.4p is uniquely solvable in the weighted Sobolev space 
Wj.Qj'(r2) as a perturbation to the background self-similar transonic shock solution. 

In particular, the linearized problem to (j3.12p . ()3.2p . and ()3.4p in the iteration 
J is 

Ao d^6U^ + Ba dr/SU^ + C{t]) SU^ = F{6U; ^) in Q, (5.3) 

Sv^ — ip'{x{^,(_tanLUo)) Sui. ~ go{5U;ip) on Fq, (5.4) 

aSu^, + (3Sv^, ^ gi{SU;il!) on Fi, (5.5) 

where Aq — A{Uq) and Bq = B{Uo) for the background solution Uq = Uq{6) = 
{uo,vo){0) between Fq and Fi described in Section 2, and 

a = — ([/o(c.i);l,0), /3=_([/o(c,i);l,0), (5.6) 

for 

G (f/; [/-):= HM + HM- (5.7) 
We denote this linearized problem as a linear operator ^ : 6U {F;go,gi) for 
F = {hj2). 

Since {du,dv)p = —p^^'^[u,v), we have 

duG (U; U^) = [pu] + {p + udup)[u] + vdvp[v\ 

= [pu] + {p- u^p^-'<)[u] - uvp^-^[vl 

d,G {U; U-) = ud,p[u] + [pTi\ + {p + vd^p)[v] 

— —uvp'^~'~'[u] + [pv] -\- {p — p'^~'~')[v\. 



as 1/ ^ 0, 



Then 

a^p{u-^ + {u- 1) (1 - u^p^-'') - uv'^p^--'^ = 0(1) 

(3 = p{-{u- l)uvp^-^ + v + v{l~- v'^p^-'^) ) = 0{l)v^ as 1/ ^ 0, 
where 0(1) depends only on 7 and h. Then 

\-\ = 0{l)u^ asi/^0. (5.8) 

a 

Therefore, there exist constants j/q and Eq such that, for any Q < v < and 
< e < eo, 

||(Ao -A,Bo- B)\\c„(i-i) + \W - tantJollcj^^^(ro) + |^ - cott^ij < e. 
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where A and B are the matrices in ()4.16p and e is the constant in Proposition [4?3l 
li F e iW°{^(n))^ and € wl~^^^''''^{Tj),j = 0,1, by Proposition SSI there 
exists a unique sohition (5?7» e (W^qI (fl))'^ to the Unearized boundary value problem 
dSS])— JESl) such that 

\\SU4wl^,^n) < K ( ll^^llw«-(si) + E llff^lU-^/'-cr,) ) ' (5-9) 

where K is independent of (5U^, F, go, gi), but depends only on b and 7. 

With the linearized problem, we will start the iteration scheme with F = (/i , /2)^ 
and gj,j = 0, 1, that take the following form: 

F{SU; V) :=C'(r/; ~ A{U; V-) £'2?/ - B{U; i') DiU 
+ Ao d^SU + Bo drjSU + C{tj) SU 

- {A{U) d^U + B{U) d,jU + C{r]) U) , (5.10) 
go{SU;ip) ■.=uo{(p'{x{£,,^ta.nujo)) - tantJo), 
gi{SU; V) + (3 5v - G{U] U'), 

where 11 — 110 + . For simplicity, write Sip = ip — cotwi, Stp' — tp' — tantJo, and 

/l = ^-T « + 1 - —)D2U - —D2V 



5V 



; - — 1 - -^)Div 



1 + tan Wo 

+ (1 - 4)9€<5u - ^ (9,<5u + d^bv) + (1 - %)dr,bv + - 

Cq Cq Cq 77 

(111 WV 1) u \ 

(1 - - ^ (^"j" + ^««) + (1 - ^)^';^ + -) ' 

72 ~ rD2V rDlU. 

1 + tan LUoSip 1 + tan o^o (5'0 

6. Proof of Main Theorem I: Fixed Point of the Iteration Map ^ 

In this section, we first prove that there exists a unique fixed point of the iteration 
mapping ^ introduced in Section 5. To achieve this, we prove that ^ is a well- 
defined, contraction mapping. 

We will need the following lemma. 
Lemma 6.1. Suppose that h{0) ~ and |l^'(e*)||^g(][j-) < 00. Then 

II^IU,(«)<i^||/.'(e*)|L,(„), 
that is, \\h\\-^o,q jjj^^ < K \\h'\\^^o,qf^^^y Moreover, for any constant s 7^ 0, we have 

\\H^^ )|L<!(R) = )|L9(R) ■ 
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These can be seen by the fonowing direct calculations: 



J-oo 6 Jo X Jq X Jq 

pOO 1 pi pi pOO 1 pi pOO 1 

</ -/ \h\sx)\''dsdx^ / -\h'{sx)\''dxds^ / -\h'{x)\''dxds 
Jo ^ Jq Jo Jo ^ Jo Jo ^ 



-\h'{sx)\'^dx= I |/i'(e*)|''dt, 

X 



and, for any constant s ^ 0, 

/OO /"OO /'OO -1 

|/i(se*)|''di= / |/i(sa;)|''d(lna;) = / — \h{sx)\'^ d{sx) 
OO Jo "'0 ■^^^ 

/•OO 

|/i(y)rd(lny)= / |/j(e*)rdt. 



6.1. Well-definedness of the iteration mapping ^ . We first show that there 
exist positive constants M and Eq such that, for any < £ < Sq, ^ is well-defined 
in Omc with the help of estimate ()5.9p . 
By Lemma l6.H we have 
II 

= 11^0 -^{ip{x) - a; tan Wo + tan wo(V'(?7) - ^cotwi)) ||^^o,,(.j^^ 

^_ /|i a;(v3(x) — ajtanwo) II n tanwo('0('7) ~ ^cotwi) , 



= 0{l)K{u;o)i^T^ ( II VIU°„«(R+) + ll'^^llvv,%''(K+)) 
= 0{l)K{ujo){l + Ms)!^^ e. 
Similarly, we have 

II II ^V-' 1 / \ o II 

^DiMo wo.-Jroi = -^-{~ sin9 + ta,nijocos9)deuo\\ o.g,ry) 

' II "ni'."'' II 1 ^ f H-n Ai/i r ii"nii."^ 



l + ta.nLUo6ij " 1 + t&nujo Sip r 



.D,{6u)\\^o.,,„, <KMsMe\ 



"l + tanc.0^^ ^llH'aTlf^) 
II (1 - )^?^ - (1 - 7:j)^«"IL°f (a) = II (72 ~ 7!) (^C"o + d^Su) llwo-cn) 
< KMv^e. 

The other terms in the expression of F can be estimated analogously. Hence, we 
have 

1 1 1 1 vyo. (H) < [Msv^^ + Mv^ + MMse^. (6.1) 
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It is easy to see that 

||.9olU^i-i/,.,(r„) < Ku^ ||VliH.i^-i/.^,(r„) < Kv^^e. (6.2) 
Furthermore, we have 

||gi|!tyi-i/9.<!(ro) - ^h^Wwl^^in) ' (6-3) 

and 

gi = {G{U;U^)~G{U;U-)) + {aSu + pSv - G{U;U^)). 
With a direct calculation, we have 

{du-,dv-)p~ = -{p~f~''{u~:V~), 

(5^ , e )p- = -{p-?-H^ + (2 - 7)(P-)'-^ 1 + (2 - 7)(p-)'-^(«-)'), 
Then 

9„-G = [u]u-v-{p-r-' - H - H(p" - (v-np-r-^), 

and 

a^G - - (Mu- + Ht;-) dl-p- +2{p-- {u-f{p-f-'<) + 2[u\u'{p~f-\ 
dl-G^ - {[u\u- + [v]v-) dl-p- +2{p-- {v-f{p-?-^) + 2[v]v-{p-f-\ 
du-v-G= - {[u]u- + [v]v~) d^-,,-p' - 2u-v-{p-f-' + {[u]u- + [v]v-) {p-f~\ 

a„„-G = -p-p- + u^p^-^ + (u-f{p-f-\ 

a^„-G = uvp^-'^ +u'v-{p-f''^ , 
5„„-G = + u-v~{p'f~'^, 

a,„-G = -p - p- + ^^2-^ + {y-f{p-f-^. 

1 

Since £ <C v^-^ , wc obtain that, for any 1 < 7 < 2, 

||G(C/; C/q-) - G(C/; =|| V^-G(C/; C/7)ds • 'JC/" 

<A'£(l + T-1 + 7-1 Me) < Ke, 

where = sU^ + (1 — s)Uq and K is independent of v and e. 
Analogous calculations for ViyG, ^^G, and V^G yield 

\\aSu + /3Sv~ G{U- U^)\\^.^^^-^ = \\l J^' %G{Us;U^)6Uds SU\\^^.^^^^ 

where Us ~ sU + {1 — s)Uo and K is independent of ly and e. 
Hence, by (|6.3p . wc have 

||gi(5(7;^0IUi-i/...(r,) + (6-4) 
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where K is independent of v and e. 

Therefore, we can choose > Q and eo > sufficiently small such that, for any 
Q < V < and < e < eo, 

¥U4^i^f^^)<ke, (6.5) 

where K is independent of e, and but depends on 7 and uj{){h). 
Hereafter, we fix M = K . Then the mapping ^ is well-defined in Omc- 

6.2. Contraction of the iteration mapping ^. We now show that, for 

5Ui,j = 1, 2, we can choose i^o > and Eq > sufficiently small such that, for any 

< < 1^0 and < e < Eo, 

W^isu') - < \ W - ^V'Wwl^'^m ■ (6-6) 

Noticing that ,7{SUi) = {F; gi, g2){SU^ = 1,2, wc have 

< /V ( \\F{SU') - F{SU')\\^o.,^^^ + J2 hASU') - 5,('5f/^)|Ui-v,.,(p^.) 

(6.7) 

where K is independent of v, cr, (5f7^, and 5Ui,j = 1, 2, but depends only on iVo(&) 
and 7. 
Since 

II " ^) - ^"'^ ll^.f (^) - 11'^"' ~ ^"'llH^<;,f (O) ' 

" -Di{Su' - <5ui)|l „,,,„, < KMss \\6u^ - Su^l, 



" l + tano;o5^ " ^^^'^^^ lll^(o,'(f2) ' 

< \\{Ao - AiU'))ddSU^ - SU')\\^o^,^^^ + \\{AiU') ~ A([/i))95C/i||^o,,(^) 

< K{Me + 0{l)i^^) \\SU' - ^t/i^^i^,,(n) ' 

and analogous estimates for the other terms of F{5U^; tp) — F{SU^; ip), we have 

\\F{6U^) - F{SU')\\^^^,^^^^ < Ke{0{iy^ + Me + Mge) \\5U^ - 5U^\\^,^t^^-^ ■ 

™(6.8) 

Obviously. 

go(<5f/';V')-5o('5f/^V) =0. (6.9) 
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Moreover, 

= a{5u^ - 6u^) + /3((5i;2 - 6v^) - {G{U^] U-) - G{U^;U')) 
= {a{6u^ - Su^) + P{Sv^ - Sv^) - {G{U^; U^) - G{U^;Ua)) ) 
+ ( (G(C/2; U,) - G{U^-U^)) - (G(C/^ U-) - G{U'-U-))). 

Then, an analogous calculation for gi above in verifying that is a well-defined 
mapping in Omc yields that 



|52(C/';V^)-52(C/^V^)| 



w, 



< KMe WSU"^ ~ <5t/M 



(6.10) 



.(0) - II IllVjofm ■ 

Choose I'o > and > sufficiently small. Then, for any < e < Eo a-nd 
< 1/ < t'oi estimates (|6.7p - (|6.10p imply that (|6.6p holds, that is, ^ is a contraction 
mapping in Oms- 

7. Proof of Main Theorem II: Fixed Point of the Iteration Map 

In this section, we prove that there exists a unique fixed point of the itera- 
tion mapping introduced in Section 5 by showing that is a well-defined, 
contraction mapping, which completes the proof of the main theorem. 

7.1. Well-definedness of the iteration mapping ^s- Let ^s{ip) ~ V'*- Write 

[v]{l — tan Wo cot wi ) 



Then 



[v\{l — tan cot wi ) 



u] + Lp'{x)[v\ 

[v\{\ — tan Wo cot loi) 



a„* = - 



[u]{l — tan Wo cot uji) 



{\u]+^'{x)[v]f ' 
Thus, by (|5.2|) . we obtain 

t/i* — cot LOl 



(M + ip'{x)[v]Y 

[u] ( 1 — tan Wo cot wi ) 
{[u]+v'{x)[v]f 



< KMe, 



(7.1) 



where if is a constant independent of v and e. 

We choose Ms = KM hereafter. Then is well-defined in Smse in the case 
that the positive constants v and e are sufficiently small. To complete the proof, it 
suffices to verify that is a contraction mapping in Y,Mse- 

7.2. Contraction of the iteration mapping ^s- Let ^s(V'"') = — lj2. 

Then we have 

r ^(SUj) = iF;go,gi)iSU,;^n. J = 1,2, 
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Thus, we obtain 

||(5[/2-5C/i||^yi,,(f,) 

< /V ( \\F{SU2;i'') - F{SU,;ij')\\^o.^^^^ 2) 

where K is independent of SUj and ^p-' ,j = 1, 2, but depends only on LOo{b) and 7. 
Since 



and 



■~y{r];ip^) 7/(7?;?/'^)^ W yimip^)^ 
1 1 tanwo(7/'2 - -01) 



we have 

II JilU«-(n) < + Me) Wi^^ - ^Ih/o,,^^^) + KMse \\Sv2 - Svi\\^i^,^^^ . 

Set 

J2 '■= —D1U2 H —DiUi 

1 + tan luq S^P^ 1 + tan luq Sip^ 

— : — D1U2 H — —Di (5ui — 5U2) ■ 

1 + tanwo^V'^ l + tanwoi^V'^ 1 + tan '^t/'^ 

An analogous calculation yields 

ll-^2|Uo.,(n) <K{v^^ + Me)||7^2 _ ^ j^^j^^ jj^^^ _ SuiW^i^^^^^^ . 

Set 

J3 :=(1 - 4)d^{5u2 5u,) ((1 - 4)d^U2 - (1 - 

'-I ^0 ^2 ^1 

Then, as the calculation for ^ , we have 

II J3|lvy«-<!(j2) < if (i^^ + Me) \\Su2 - SuiW^i^,^^^ . 

Analogous calculation for the other terms of F{5U2;i^'^) — F{5U\\^^) finally leads 

to 

||F(5t/2;V^2)-i^('^t/i;7/>i)||^o,,(^) 

< K{v^^ + A/e)||7^2 - 7/,i||r^ + K{Mse + Me + t^^) \\5U2 - 6Ui\\^i^,^^^ . 

(7.3) 
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Since 

=uoip" ('(/'^(Ctanwo) - ^'^(Ctanwo)), 

and 

■^ip'{x{^, ^ tan Wo; V'"')) = y^"{x) (^d^x + ta,nLUodrfx) ~ ip" {x) (l + tanwo^'f/'"') , 
we have 

||go(<5t/2;V')-5o(<5C/i;V'')|Li„-i/,..(r„) <^i^^^llV''-^'l|ri, (7.4) 

where K depends only on wo(6) and 7, but independent of -0-', z/q, and Eq- 
Furthermore, we have 

5i(<^C/2;V^')-.9i(^C/i;V^') 

^a{5u2 - <5^ii) + /3(5w2 - - (G(C/2; [/-(ri;^'^)) _ q^^^. ^-(r,;^!))) 
=a(5u2 - 5ui) + /3(,5v2 - - (G([/2; [/-(ri;^^)) _ ^-(p,; v?))) 

+ (G(;7i; ;7-(ri; V')) - G([/i; [/-(Fi; 0i))), 

where [/^(Fi;-!/;-') [/~ (77 cot oji, 77; ■;/'-'), J = 1,2. Notice that 
G((7i;(7-(Fi;V/))-G(C/i;C/-(Fi;7A')) 



(a„-G([/i;C/7) (5u-(Fi;7/>^) - ,5u-(Fi; tA')) 
+ 5,-G([/i;[/7) (5z;-(Fi; V') -<5z;~(Fi; V')) )ds, 
where [/- = s?7-(Fi;V'^) + (1 - s)C/- (Fi; V-i), and 

^C/-(Fi; 7^2) - ^t/-(Fi; V') = Jf/-(V''(7?), ry) - 5U-{i^\ri), ri) 

d^U~{stl;^ + (1 - s)V'\ ry)ds (i/;^ ~ . 



Then an analogous calculation as for ^ yields 

< KiMe \\SU2 - SUiW^i^,^,-,-^ + K2e\\i;^ - V^^liri, 

where Ki and K2 depend on tooib) and 7. 
Then, by we have 

ll'5^2-^C/i|lvK^i^.«(a) 

< K{i^T^ + Me)\\i;^ - i;^\r, + K{Mse + Me + f^) \\SU2 - 6Ul\\y^,l, 



(7.5) 



(0) • 

(7.6) 
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Choose and eo sufficiently small. Then, for any < i/ < i/q and < e < £o, 
we have 

\\SU2-SU,\\^u.^,,^<K{,.^ +Me)U'-i,^\\r,. (7.7) 

Thus, 

< K \\SU2 - SUiWy^.u.^,,^ + K \\6U-{Ti;ip^) ~ dU-iTi;^p')\\^^ 

< K{iyT^ +Me + s)\\iij^ -^'^IIti 

where, for the last inequality, we have again chosen vq and £o to be sufficiently 
small. 

This implies that is a contraction mapping so that it has a unique fix point 
in YiMsei which completes the proof of Theorem 13. II 

Appendix: A Fredholm-type Theorem 

To be self-contained, in this appendix, we give a proof for a Frcdholm-type the- 
orem. Theorem A.l, a special case of Theorem 4.1 in Maz'ya-Plamencvskii [28] . 
following their ideas. Consider the boundary value problem of an elliptic equa- 
tion of second-order in an infinite strip Q := {(i, x) : x € / := (.xo, xi),t G M} with 



boundaries Eq = = ^o} and Si = {x = xi}: 

Lip := dtt^ + dxxtp + dtip + a{x)dx^ ^ f inQ, (A.l) 

Bqip dxip = 50 on Sq, (A. 2) 

Boip := dxip ^ gi on Si, (A. 3) 



where a{x) £ C\l), f e W°:_\^{g), and £ wl_^^^'^''^{R),j = 0,1. We assume 
q > 2 since only this case is really used in this paper. Obviously, the operator 
(L; Bo,Bi) of the boundary value problem (|A.ip - (jA.3|) acts continuously from the 

space w^:\^ig) to w°:_\^{g) X (Ty(i_-y«'«(M))'. 

Consider the boundary value problem with a complex parameter A on the interval 

/: 

d^^ip + a{x)d^ip+ {-X^ +iX)tp ^ f in/, (A.4) 

dx^ = go x = xo, (A. 5) 

dxf = gi x = xi. (A. 6) 

For all A, with the exception of certain isolated points, (|A.4p - (jA.6|) has a unique 
solution If £ W^ P. The exception isolated points of A are called spectrum of 
problem ((X4)) ~ (|X6)) . 
Then we have 
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Theorem A.l. // the line IniA = /? does not contain the eigenvalues of prob- 
lem (|A.4|) - (jA.6p . then the operator {L; Bq, Bi) of problem (jA.ip - ()A.3P realizes an 
isomorphism: 

Moreover, the solution ip € W'^'^^ifj^ of (jA.ip - (|A.3[) satisfies the estimate: 

Remark A.l. In the case p = 2, this assertion is well-known (cf. [29]). In this case, 
a solution in the class W'^^^^IG) can be represented in the form 

^{t) = I e-'^'R{\),^t^x {/; 30, 9i} dA, (A.8) 

V 27r JlmA=^ 

where i?(A) denotes the inverse operator of problem (|A.4[) - (|A.6p and J^t-tX is the 
Fourier transform with respect to the ^-variable into the A-variable. If it is addi- 

Q2 ]^ 1/22 

tionally assumed that / G ^(^fj-^^iG)} 9j G ^(/3i) ' (^) ^^^'^ that, in the closed strip 
between the lines Im X = (3 and Im A = /3i , there are no points of the spectrum of 
(fO)) - ([X6l) . then the function ip defined by (fO|) belongs to w'^l^^^{Q), and 

To prove Theorem lA.il we need two lemmas, which are all in [28]. 

Let i2^o, ■^1; and ^2 be Banach spaces of functions on M, in each of which 
multiplication by scalar functions in C^(M) is defined. Let {C,k}°^oo ^ partition 
of unity on R subordinate to the covering of M by the intervals {k~\)5 < t < (fc+1)^, 
where 5 is a fixed positive number and Cfe & C°°(R). Suppose that the norms ||-|j^ 
in the spaces = 0, 1, 2, possess the following properties: For p G [1, 00], 

00 ^ , 

Ci\\ul<[ IKfe^llo) '<^^2||^.|lo, (A.IO) 

A;— — 00 
00 ^ / 

\\vh>c[ \\<^<) (A.ll) 

A;— — 00 
k— — oo 

Lemma A.l. Let : J^i — > ^0 be a linear operator defined on the functions with 
compact support and such that, for some e > and any integers m and k, 

\\Qk^{Qmv)\\^<Cc-\-^-''\'\\Uv\\^ foranyvGJ^i. (A.13) 

Then 

(i) For all v G jz/i with compact support, 

\\^v\\,<C\\v\\,, (A14) 
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where the constant C does not depend on v. 

(ii) Let ^2 C ^0 • Suppose further that, for all functions v in with compact 
support on M, 

\\Ck^v\\^<C{\\auv\\^ + \\ak:'3^v\\,), (A.15) 

where au = Cfc-i + Cfc + Cfc+i; k = Q, ±1, Then 

\\^v\\^<C\\v\\,. (A.16) 
Proof. According to (jA.lOp and (jA.lSp . wc have 

i/p 

k— — OQ m— — oo 
-, / oo oo T / 



io = lh^( E Mllo<^^( E II E a^(c™^')||^) 

?7i — — oo /c — — OO rri — — OO 

oo oo -. / oo oo 

<^( E ( E W^'^nc^mi) <c{ E ( E k^vw^y 



fc— — oo m— — oo A:— — oo m— — oo 



Since the operator of discrete convohition with kernel {e ''^ £^cts continuously 

in Ip, it follows that 

i/p 



m— — oo 



The last inequality, together with ()A.lip . leads to (|A.14p . 
Furthermore, by (|A.12|) and (|A.15|) . 

ii^«ii2<c(EiiCfc^Hi^)'^'<t^(Eii'^^^c)'^' + ^^(Eii'^^^«iio''^' 

Using the definition of tTfc, (|A.10[) . and (|A.11[) . we obtain 

i/p 



Then we apply (|A.14[) to arrive at the result. □ 

Lemma A. 2. Suppose the supports of the functions f and gj,j = 0,1, are con- 
centrated on the set {{t, x)(lQ:m^l<t<m + l} (m an integer), and f G 
W°'P{g) n W^-^iG), gj e W^-^/P^P{R) n W^-^/^^^{R), j = 0,1, for p > 2. If 
the line Im X = f3 does not contains the eigenvalues of problem (|A.4[) - (jA.6|) . then 
the solution ip e W'^^'^{Q) of problem (jA.ip - (|A.3P satisfies the estimate 



e''*O^IL.(,)<Ce"l™-'l^ 11/11 5: , (A.17) 



.=0,1 



where e is a positive number and {Ol'^^oo partition of unity on M subordinate 
to the covering of R by the intervals I— l<t<l + l. 
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Proof. Denote by Mp.p the term in the parenthesis of the right-hand side of (jA.17p . 
Using ()A.9P for any /3i g (/? — e, /? + e), we obtain 

( h'Mll^n'^ty^' <Ce('3-^^)'M2,ft < Ce(^-'^^)'Mp,^, 

VJ;_2 ^ ' J (A.18) 

<Ce('3-/'i)('-™)Mp,0. 

Applying results in [1], we find that the solution ip e W^^^iG) of (lAl|) - f01) has 
the following local estimate: 

\\vM\w^.P(g) <c[\\V2f\\w«.p{G)+ \\^^Sj\\w^-^/p,pm + \\^M\L^(g)), (A.19) 

j=0.1 

where rjsit) = vit/s), rj G C^{-1, 1), and 7](t) = 1 for \t\ < 1/2. 

In the case |m — ^| < 2, the local estimate (|A.19|) leads to the estimate 



O^LHG) ^ Cc^^-^^^^'-"'^M,^p + c[ ^^'||e^V|r^.(,)dt)'^'. (A.20) 



If |to — ^1 > 2. then, by (|A.19p . the last inequality remains valid even without 
the first term on the right-hand side. Combining (lA.lSp with (jA.20p . we obtain 

VCM\lhq) ^ C'c('3-^^)('-")Mp,0. (A.21) 

Setting /3i = /3 + £ for m < I and (ii — (3 ~ e for m > I, wc arrive at (jA.17p . □ 

Now we prove Theorem lA.il 

Proof of Theorem lA.li Existence. It suffices to prove (|A.7p for a solution ip G 
n W^if^iG). Let ^ be the inverse operator of problem (|XT|) - ((X3l) defined 
by ^KM on the space W^^^^iG) x {W^^^^'^ (R))^ . We set 



Lp(g) ' 



\\{f;9o,9i}\\^,--\\f\\w°^''{g)+ J2 ll-9jll<,-:i/'''''(K)' 



.=0.1 



Lemma IA.2I and (jA.19[) ensure that the hypotheses of Lemma lA.ll arc satisfied. 
Therefore, for the solution (f = ^ {f; go, gi}, we have (|A.16p or, equivalently, 

(TOl) . 

Uniqueness. Let ip he a solution of the homogeneous problem (jA.ip - (jA.3p in 
W^f^iG). We set Gs := {{x, t) e G : s < \t\ < s + 1}, s ^ 1,2,. . . , and introduce the 
sequence of functions tps G C^(R), ips — ^ for |<| < s, tps = for |t| > 1, and 



< C, j ~ 1, 2, for some constant C < oo. Then the function tpsf satisfies 



( !0|) - ((A31) . where / and gj are functions concentrated in Gs ■ Since p> is a, solution 
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of the homogeneous problem, by the local estimate (|A.19p . we have ip £ W^ioc (^) 
and 



Since the right side of this inequality tends to zero as s oo, it follows that ip ~ 0. 
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